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Based on an idea due to Bardeen and Pearson, we formulate the light-front Hamil- 
^ . Ionian problem for SU (N) Yang-Mills theory in (2 + l)-dimensions using two continuous 

space-time dimensions with the remaining space dimension discretized on a lattice. We 
^ ; employ analytic and numerical methods to investigate the string tension and the glue- 

;h ' ball spectrum in the N ^ oo limit. Our preliminary results show qualitative agreement 

^ with recent Euclidean lattice Monte Carlo simulations. In the following, we attempt to 

give a more pedagogical introduction to the idea of the transverse lattice; more detail 
may be found in Ref. |jl||. 



MOTIVATION 

There exists a large gap between the quantum field theory of QCD with its many 
successes in the context of perturbation theory and experimental observables associated 
with QCD bound states: the hadron mass spectrum, structure functions, form factors, 
et cetera. There has been little progress during the last 20 years in building a bridge 
between the field theory and these non-perturbative properties of the hadron spectrum. 

* Presented by B. van de Sande 
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Let us contrast two attempts to bridge this gap: the Euchdean lattice Monte Carlo 
(ELMC) approach and light-front field theory where one uses Hamiltonian techniqes 
on a theory which is quantized on a surface of constant = + x^)/\/2. As far as 
progress is concerned, ELMC is much further along, benefiting from a large research 
effort since the mid 1970's. In contrast, the light-front approach is not so far along; 
most research effort has occured since late 1980's. For ELMC, further progress is limited 
mainly by the speed of available computers. For light-front, progress is currently limited 
by conceptual issues, exempli grati renormalization. Even when successful, ELMC is 
not able to measure many interesting observables, such as structure functions, directly. 
For light-front field theory, physically interesting observables are quite easily calculated 
from the bound state wavefunctions. The approach that we will discuss here, the 
transverse lattice, uses ideas from both the light-front and lattice approaches. 

Instead of solving the full theory of QCD, we will focus on a particular model: 
SU {N) Yang Mills theory in 2+1 dimensions in the N ^ oo limit. Why 2+1 dimensions 
instead of 3+1 dimensions? Aside from the obvious fact that it has fewer degrees 
of freedom, one should note that the theory is super-renormalizable. In the context 
of lattice calculations, this means that there is no critical point and the associated 
'critical slowing' is absent. Consequently, excellent lattice spectra are available for 2+1 
dimensions which one can use as a comparison. Why large A^? The N ^ oo limit 
allows considerable simplification of our computational problem; in addition, the lattice 
data indicates that 1/A^ corrections to the low energy spectrum are quite small. 



THE TRANSVERSE LATTICE 

A number of years ago Bardeen and Pearson ^ formulated a light-front Hamil- 
tonian version of lattice gauge theory, which makes use of the fact that two components 
of the gauge field are unphysical. In this approach two spacetime dimensions are kept 
continuous = (x° ± x^) / while the remaining transverse spatial dimension x^ 
is discretized on a lattice. Lattice sites are labeled by integer-valued index i and the 
lattice constant is a. Following Bardeen and Pearson, we associate longitudinal gauge 
fields Af(x^,x^) with lattice site i and link variable Ui{x~^,x~) with the link between 
sites i and i + 

X X . (1) 

i i + 1 

Note that Af and Ui are N x N color matrices where Af is Hermitian and traceless 
and Ui G SU{N). We can define the color Maxwell tensor 
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where a, (3 G {+, — } along with a covariant derivative 

D^U, = (9" + lAU) U, - tU,A1^, . (3) 
Using this notation, one can write down an action, 

A = aJ2j dx+dx- Tr|-^F„^,,i^"^ + ^^D^U^D'^U}^ , (4) 

which has several important properties: 

• If we take the naive continuum limit a ^ where Ui = exp{—iaA}), we recover 
the continuum Yang-Mills action. 

• By construction, this action is 100% gauge invariant. 

• It has an automatic confinement mechanism. We will say more about this later. 
Linearization 

In an ideal world, Eqn. ^ is the action one would quantize. However, we are 
faced with a problem. The field theory on each link of this transverse lattice is an 
SU{N) X SU{N) non-linear a-model. Upon quantizing this theory, we must enforce 
iV^-nonlinear constraints on each link. Although some attempt has been made to do this 
quantization success has remained elusive. Instead, we perform the linearization 
suggested by the original authors |^. In this scheme one replaces the unitary link 
variables Ui with N x N complex matrices^ Mj, that is, Ui — > \/2ag^Mi. 

The obvious next step would be to include an effective potential to enforce the 
constraint \p2ag^Mi G SU{N) dynamically. For instance the termQ 

A,Tr|(M;M,-2a/)'| (5) 
= XcTr[MlMiM}Mi} - Aag^ Tr {M^Mi} + (2a/) ^ (6) 

is minimized precisely when \f2ag^Mi G U{N). One could imagine adding such a term 
to the Hamiltonian, taking the Xc —>■ oo limit, and recovering the transverse lattice 
action (^). Unfortunately, a closer inspection reveals several problems: 

• The second term in is a mass term with negative a (mass)^ coefficient. Quan- 
tization of the theory breaks down for negative (mass)^ for the same reason that 
it breaks down in the 't Hooft model and the spectrum becomes unbounded 
from below. Presumably this negative (mass)^ term is a signal for the presence of 

*The 0{n) cr-models are a good example of where this hnearization is known to work. 
^For large N, we can ignore the distinction between U{N) and SU{N). 
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spontaneous symmetry breaking. One could imagine solving the broken phase of 
the theory using zero mode techniques 0. The resulting 'shifted' theory would 
have some complicated effective potential representing the effects of the sponta- 
neous symmetry breaking; this effective potential is undoubtedly different than 
the above form (|^). 

• Tadpole contractions associated with four-point interactions, such as the first term 
in d^), produce a divergent shift in the mass term Tr{M/Mj}. Renormalization 
leaves the mass as a free parameter in the Hamiltonian. 

• A successful renormalization group analysis would introduce couplings between 
neighboring links. For instance, one expects terms associated with the constraint 
2ag^MiMi+i e SU{N). 

• From a more practical viewpoint, implementing the Xc ^ oo limit forces half of 
the dynamical degrees of freedom at each link to decouple from the theory. This 
is not good news for numerical calculations where the computational difficulty 
depends critically on the number of dynamical degrees of freedom. 

As a consequence of these considerations, a first-principles construction of the correct 
effective potential Vi is not a simple matter. Instead, we take a more pedestrian ap- 
proach and include in Vi all operators up to fourth order in Mj containing one color 
trace 

Vi = fi''Tr{M,Ml} + ^TT{MiM}MiM}} 

+ ^Tr{M,Mi+iMl,Mj} (7) 

and try to determine the associated coupling constants empirically. Note that the Ai 
term is local to one link and the A2 term acts on two adjacent links. 



Quantization 



Next, we quantize the theory, write down the Hamiltonian P~, and construct a 
basis of states. Further details may be found in Ref. At each lattice site i we have 
a 1-1-1 dimensional gauge theory with conserved current 

= ilM.B' M} + mI^B" Mi_^ . (8) 



We set d^Af = by choice of gauge and throw away the associated dynamical zero 
mode / dx~ Af . The A~ field obeys the equation of motion 

{d-fA- = ^^Jt. (9) 
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As evidenced by the absence of time derivatives in Eqn. (|]), A~ is not a true 
dynamical degree of freedom but is constrained. We solve the constraint equation (||) 
for A~ and remove it from the theory. At this point, quantization of the theory is 
straightforward. The momentum conjugate to Mi{x~) is d-Mj{x~) and we impose the 
usual equal commutation relation^ 

'M,{x-),d.M]{y-)] = '-S,,S{x- -y-) . (10) 
The longitudinal momentum operator 

= 2j2j dx-Ti[d-Mid.Mj} (11) 

i 

and Hamiltonian 

p--Eid.-(-f^Tr[jtji^jt}^v:^ (12) 

generate translations in the x~ and x~^ directions, respectively. 

Finally, we construct a basis of states. The zero mode of the A~ constraint equation 
(^ generates the Gaufi' law constraint 

= [ dx-Jt (13) 



which we must impose on the basis of physical states. That is, physical states must be 
color singlets at each lattice site. We construct a basis of closed color loops: 



Tr{MiM/} |0) t % 

i:i[MiM}MiMl] |0) 



X 



Tr{M,M,+iMiiM;} |0) > ^ x ^ (14) 



Tr{M,M,+iMi.iM,+iMiiM;} |0) 



3< 



et cetera. i i + 1 i + 2 

where the x~ co-ordinate of each link field remains arbitrary and the number of links 
Mj must equal the number of anti-links Mj. Since the loop- loop coupling constant is 
non-leading in A^, we do not include states with more than one color trace in the Hilbert 
space; we deal with a free string theory.[] Incidentally, if we do introduce a state that 

*The factor of 1/2 comes from the Dirac procedure for constrained systems. Also, we drop the 
constrained zero mode Jdx~Mi. 

*In fact it is possible to choose Vi such that one has a theory of free bosonic strings ||^. 
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is not a color singlet at some site, exempli grati 

Tr{---M,M,+2---}|0), (15) 
we find that its energy diverges. 
Confinement 

This theory has built-in linear confinement. Consider two test charges separated 
in the x~ direction. The first term in Eqn. ( [T^ ) acts as a linear confining potential.^ 
Now consider two charges at lattice sites i and j. Due to the GauB' law constraint (IT^) 
we must construct a color string of at least link fields between the two charges. If 

/i^ > 0, there is some minimum energy associated with each link field and the transverse 
string tension is nonzero. As we shall see in the next section the resulting potential is, 
in fact, linear. Similar arguments hold for the original action (^. 

Numerical Techniques 

We solve the spectrum on the computer using DLCQ techniques [^]. For the x~ 
coordinate, we impose anti-periodic boundary conditions Mi{x^) = —Mi{x~ + L) and 
use a momentum space representation. For integer valued cut-off K = LP^ /{2'n'), 
momenta are labeled by odd half integers Km G {1/2, 3/2, . . .} where J2m = K. This 
yields a finite basis of states. Also, we employ a truncation in particle number. 

The first term of the Hamiltonian (0) dominates the behavior of the theory. Thus 
it is natural to choose units such that the associated coupling constant g'^N/a is set 
equal to 1. This choice of units is assumed in the following. 



STRING TENSION 

Next, we introduce a method for measuring the string tension in the x^ direction. 
Consider a lattice of n transverse links and periodic boundary conditions. We construct 
a basis of states that wind once around this lattice and calculate the lowest eigenvalue 
of the invariant (mass) 2 operator M^]^) = 2P+P-|^). The cont inuum limit string 
tension is 

lim — - — , na fixed. (16) 
a An 

This definition is be equivalent to the standard definition of string tension using two 
test charges if the charges are placed sufficiently far apart. If we plot vs. n in 
Figure we obtain a quadratic corresponding to a linear confining potential. 

In fact, we find that the string tension decreases with decreasing /x^, leading us to 

*Note that (5_)-2/(cc-) ^ S dy-\x- - y-\f{y-)/2 in Eqn. 
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1 2 3 4 5 



n 

Figure 1. Lowest eigenvalue vs n for states that wind once around 
the lattice. Here, K = 10.5 or 11 {{n + 4)-particle truncation), ij? = 0.1, 
Ai/a = 1, and X2/a ~ —1. Also shown is a fit to a quadratic. 

the conclusion that /x^ is a measure of the physical lattice spacing. The continuum limit 
is partially fixed by requiring vanishing A(M)/A(n); see Fig. |[ Above the surface, the 
string tension is positive corresponding to nonzero lattice spacing. 

In addition we have an analytic ansatz for these states that wind once around the 
lattice that is valid in the X2/CL < region of parameter space. It agrees well with the 
numerical results. 



SPECTRUM 



Before we look at the details of the spectrum, we can make some statements about 
the allowed region in parameter space. A numerical estimate of the 'edge' of this region 
is plotted in Fig. Comparing Figs. ^ and H, we see that a continuum limit with a 
non-tachyonic spectrum occurs only for the "wedge shaped region" — Ai < A2 < Ai/2. 
In this region, the continuum limit occurs for /i^ ^ to within numerical errors. Finite 
lattice spacing corresponds to /i^ slightly above the surface in Fig. ||. 



Symmetries 



The theory possesses several discrete symmetries. Charge conjugation induces the 
symmetry C : (Mj)^^ (M/)^,; where l,m & {1,...,A^}. Parity is the product of 
two refiections Pi : — > —x^ and V2 '■ —x^. In light-front quantization, Vi is 

an exact symmetry Vi : Mi mI^ while V2- x^ ^ x~ , is complicated. Its explicit 
operation is known only for free particles |jlO[, which we call "Hornbostel parity." The 



latter is nevertheless useful since it is often an approximate quantum number and can 
be used to estimate 7^2 • Given V2 and Vi we can determine whether spin J' is even 
or odd using the relation (— 1)"^ = ViP2- If rotational symmetry has been restored in 
the theory, states of spin 7^ should form degenerate Vi doublets \+ J) ±\ — J) . 
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Figure 2. Parameters such that the lowest Af^ eigenvalues are equal for 
n = 4 and 5 (see Fig. [^), where K = 10.5 or 11 ((n + 4)-particle truncation). 
This is an estimate of vanishing string tension. Also shown is a line such that 
the M"^ eigenvalues are approximately degenerate for n = 3, 4, and 5. 




Figure 3. Parameters such that the lowest eigenvalue is zero, K — 10 to 
14 with extrapolation using a fit to X]m=o '^mK~'^^'^ (6-particle truncation). 
Below this surface, the spectrum is tachyonic; above the surface, it is well 
behaved. 
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Figure 4. A comparison of our spectrum with SU{3) ELMC data in 
units of the physical string tension ||] for various \J'\'^^^. The parameters 



g^N/a = 3.90, ij/ = O.lMg^N/a, Ai 



0A87g^N, and A2 

,2 



l.lOSg^N were 

chosen by a best fit to the lattice data, x'^ = 40, where i^T = 10 (8-particle 
truncation) . Our error estimates are solely for the purpose of performing the 

x' fit. 



As with the lattice results, we use "spectroscopic notation" \J'\'^^^ to classify states. 
One expects the lowest two eigenstates to be approximately two-particle states 

Y^TT{Mkx')M,iy-)}\0) (17) 

i 

with the lowest state having a symmetric wavefunction corresponding to and the 
first excited state having an antisymmetric wavefunction corresponding to . Of 
course, these states also have 6-particle et cetera contributions. 



Results 



Ideally, we would like to predict the effective potential based on some connection 
to the continuum theory, restoration of rotational invariance, et cetera. However, as a 
first step, we use instead a best fif fo the ELMC results of Teper |0. 

An example spectrum is shown in Fig. ^. Similar spectra are found in other regions 
of coupling constant space above the 'wedge-shaped region' (Fig. |). We label the lowest 
2 and second states based on the expectation value of the number operator and 
determine (—1)"^ based on Hornbostel parity the exception is the IJT'I'''" sector 
where Hornbostel parity gave exactly the opposite of the desired results. Beyond this. 
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J is determined by a best fit to the lattice data. 

Since this spectrum is the result of a best fit, it is not very predictive. However, 
we can use the result to tell us about our model. At first glance, we seem to have a 
pretty good match. However, we note several problems: 

• The energy of the lowest state is too low. 

• The lowest 2"*"+ and states form a degenerate doublet if rotational symmetry 
is restored (as indeed happens for the lattice data). In our case, the splitting is 
large. This discrepancy dominates the error in our fitting procedure. 

Let us review the possible sources of error in our calculation. 

Large N. We compare N oo spectra to SU{'i) lattice results. However, based 
on lattice calculations for SU{2), SU{3), and S't/(4), corrections to the low 
energy spectrum are small 0. 

Finite K. Our discretization of the longitudinal momentum introduces some 
error. However, we have generated spectra for K = 10, 11, 12, 13, 14, extrapolated 
to large K (6-particle truncation), and compared to large extrapolated ELMC 
spectra. We saw no real improvement in our results. 

Particle Number. We also impose a truncation in the number of particles. We 
have examined spectra for 4-, 6- , and 8-particle truncations {K = 10), extrapo- 
lated to large number of particles, and compared to large N extrapolated ELMC 
spectra. We saw no real improvement in our results. 

Hamiltonian. The effective potential Vi that we chose (|^) did not contain any 
6-point or higher interactions. In addition, we did not include any operators 
containing multiple traces, for instance ( Tr{MjM/})^.p 

CONCLUSIONS 

We have investigated the transverse lattice model of Bardeen and Pearson [0] for 
(2 + l)-dimensions in the large- limit using linearized link variables and an empirical 
effective potential Vi. We identified a choice for Vi corresponding to vanishing string 
tension. The glueball spectrum in the vicinity agreed qualitatively with that coming 
from the presumably reliable Euclidean lattice Monte Carlo results 0. Most impor- 
tantly however, we did not see significant signs of rotational invariance which could lead 
one to conclude that the transverse gauge dynamics were correctly accounted for by Vi. 
We believe that our choice (|^ is probably too simple and that higher order terms are 
necessary to see improvement in our spectrum. 

*Although this term is generally not leading order in iV, it does act on the two particle subspace of 
the theory. 
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Future work includes the addition of more operators in the effective potentiaL Also, 
we can use our method of measuring string tension to measure the physical lattice 
spacing. This issue needs further investigation. Most importantly, we need a more 
concrete connection between our model and the continuum theory. This would allow 
us to better predict the correct effective potential. 
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